Optimal control of the local electromagnetic response of nanostructured materials: 

optimal detectors and quantum disguises 
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We consider the problem of optimization of an effective trapping potential in a nanostructure with 
a quasi-one-dimensional geometry. The optimization is performed to achieve certain target optical 
properties of the system. We formulate and solve the optimization problem for a nanostructure that 
serves either as a single molecule detector or as a "quantum disguise" for a single molecule. 



The rapid advance in fabrication techniques makes 
it possible to control the composition and geometry of 
nanostructured systems and materials on atomic scales 
[H . One of the most interesting and challenging problems 
in modern nanoscience is how to utilize these fabrication 
capabilities to control the interaction of electromagnetic 
fields with nanostructured media. A successful solution 
to this problem would make possible the design of, for ex- 
ample, high energy density storage applications, or multi- 
functional nanophotonic materials with predefined prop- 
erties [3]. 

There is, in general, a complex relationship between 
the effective electron trapping potential in a nanostruc- 
ture and the nanostructure's optical properties. One can 
usually solve this problem for a chosen geometry and 
structure of the system, thus for a given trapping po- 
tential. However, the inverse problem to find an optimal 
geometry of the nanostructure to attain desired optical 
properties seems much more complicated. In particu- 
lar, the matter-field interaction becomes complex in the 
quantum limit of very small nanostructures, containing 
only a few electrons. Quantum mechanical corrections 
and the non-locality of the electromagnetic response be- 
come significant for system sizes comparable to a typical 
electron wavelength. This condition may be satisfied for 
such small objects as quantum dots, or for the next gener- 
ation of microchips For example, for the semiconduc- 
tor GaAs with a doping level of 10^'^cm^'^ and effective 
electron mass m* — 0.07 x me, the Fermi wavelength 
Af ~ 20nm is on the length scale currently available in 
microchip design. Thus, the problem has both funda- 
mental and applied significance, and there is a strong 
demand for methods to efficiently explore the optimal 
design of quantum nanoscale devices. 

A systematic approach to the theoretical analysis of 
optimal design of potentials in quantum scattering in 
semiconductor nanodevices was first developed in 
Twelve years later, an optimal design problem with simi- 
lar objectives was considered in 0]. In a simple tight- 
binding model was considered in two dimensions with the 



aim of fitting a target density of states, and as a result, 
also describing the basic response properties of the sys- 
tem. Another interesting example of inverse band struc- 
ture engineering was considered for nitrogen impurities 
in GaP [3]. 

One of the potential applications of optimal design of 
nanostructures is to the engineering of efficient single 
molecule detectors. Single molecule sensing by means 
of optical excitation is receiving increasing attention [8]. 
In this work we will demonstrate that with the help of 
quantum design it is possible to control and manipulate 
induced local fields in the nanostructure on a very fine 
scale. We formulate and solve the optimal design problem 
for nanostructures, which may serve as single molecule 
detectors. We consider electrons trapped in the effective 
potential of the nanostructure. The presence of a po- 
lar molecule in close vicinity of the nanostructure may 
slightly modify the electron trapping potential, and as 
a result, its optical response. If a monochromatic elec- 
tromagnetic field of frequency co is used to excite the 
nanostructure followed by a measurement of the molecu- 
lar induced changes, then it is possible to use this nanos- 
tructure (or, better, an array of identical nanostructures) 
as a single molecule detector. The optimal design goal is 
to find a nanostructure with an effective trapping poten- 
tial that results in maximum contrast between the optical 
response either with or without a particular molecule at 
a given excitation frequency. Note, that the opposite op- 
timization problem, to find the effective trapping poten- 
tial that minimizes the contrast, can also be formulated. 
We also consider this problem, and refer to this type of 
nanostructure as providing quantum disguise. 

We assume that there are N^i electrons trapped in a 
nanostructure. The system is described by the Hamil- 
tonian i7„ = + Vtrapii^), with the eigenenergies and 
eigenf unctions Ei and ipi, respectively. Since inhomo- 
geneous systems with a few electrons can have an elec- 
tromagnetic response very different from the bulk, in 
our simulations we use the non-local density-density re- 
sponse function within the coordinate-space representa- 
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tion x(r, v' ,uj) : 



X(r,r',w) = ^ 



Ei — Ej — huj — ij 



(1) 



where f{Ei) is the Fermi fiUing factor and the small con- 
stant 7 describes level broadening, to is the frequency of 
the external electromagnetic field. The electron eigenen- 
ergies Ei and eigenfunctions ipi are obtained numerically 
using real-space discretization of the Schrodinger equa- 
tion. Using Eq.U) we calculate the induced electric field 
E.md(r, Lo) in the system within linear response theory 

The effect of the presence of a molecule is modeled by a 
potential which is added to the effective trapping poten- 
tial in the nanostructure Kno/(r) = — Aexp(— i?|r — rpp), 
where Tq is the position of the molecule. The Gaussian 
potential is chosen just for demonstration purposes. We 
also performed simulations using a r^^ potential with 
a cut-off, and obtained qualitatively similar results. A 
more realistic molecular potential could be deduced from 
ab initio simulations. 

As stated above, the presence of a molecule can modify 
the response of the nanostructure. To characterize the 
difference between state 1 without the molecule and state 
2 with the molecule attached, we introduce the contrast 
ratio between the two responses: 



where the normalization intensity Iq 
/y l|EL<i(r,t^)|pd3i.^ and Ei„rf(r,cj), E,|;^^(r,w) are 
the induced electric fields in state 1 and state 2, corre- 
spondingly. The choice of the normalization intensity Iq 
is dictated by the aim of finding a nanostructure, whose 
response in the state 2 generates a local electric field 
with the largest intensity. 

We assume that the effective trapping potential 
Vfrap(r) can be represented as a sum of Np model po- 
tentials Knod(r) located at r^, i = 1, ..,Np, which can be 
controlled during the material fabrication procedure with 
sufficient precision: 



(3) 



For simplicity we assume l^od(r) to be the Coulomb po- 
tential of a charge Q with a cutoff a: 



Vmod(r) = (3/|r| 
= Q/a, 



\r\ > a, 
r| < a. 



(4) 



The optimization algorithm searches for the optimal pa- 
rameters Ti to maximize or minimize the contrast ratio 
given by Eq.([2). The search for the optimal potential 



Vtrapi^) is performed using Brent's "principal axis" op- 
timization algorithm 10] for a scalar function of several 
variables. 

In order to reduce the computational complexity of the 
problem we perform the simulations considering a quasi- 
ID nanostructure geometry. We assume that spatial lo- 
calization of electrons in the other two dimensions (y, z) 
creates relatively large energy gaps in the sub-bands, so 
we can neglect excitation in these directions. The in- 
coming field is linearly polarized, with the polarization 
along the x axis. Such systems can be created even on 
an atomic scale using STM techniques [ll| by placing 
atoms of different species in a line. Let L be the length 
of the system. Then Eq = is the unit of energy 

determined by the geometry of the system. Here mg 
is the electron mass. We perform optimization of the 
trapping potential at the frequency of the external field 
hbj = Q.IEq. We assume there are N^i = 15 electrons 
trapped within the nanostructure. We set the tempera- 
ture to T = OK, and the level broadening 7 = 10~^i?o- 
For the model potential in Eq. (|4]) we use the cut off pa- 
rameter a = 0.05L and Q = 15e/Np, where e is the 
electron charge. For our simulations we set the molec- 
ular position at xq — L/2. We set the parameters of 
the molecular potential Vmoi{^) to A ~ 0.005e/i and 
B = 0.14L-2 

Figs. [1] (a)-(b) show the results of calculations for a 
non-optimized nanostructure. Fig.[l](a) gives the ground 
state electron density and the trapping potential with- 
out and with a molecule. Note, that in all the figures 
the trapping potential is scaled down by a factor of 10. 
Fig. \V[h) shows the induced field intensity |Ei„(ip for the 
trapping potentials given in FiglUa). Note, that the ad- 
dition of the potential Vmoiix) just slightly modifies the 
induced field intensity. 

Figs. [Hc)-(d) are the same as [Tla)-(b), but for an 
optimized potential. The optimization is performed to 
maximize the contrast ratio given by Eq.(l2]). Note, that 
in the optimized case, the difference in the induced field 
intensity with and without the molecule is so big, that 
we have to employ a logarithmic scale. The explanation 
of this effect is that the potential Vtrap + Vmoi produces 
an energy difference between the highest occupied and 
lowest unoccupied states that very well matches the fre- 
quency of the probe field hoj = 0.1£^o- 

Now let us consider the opposite problem of searching 
for the trapping potential that minimizes the contrast ra- 
tio in Eq. ([2]) . During the search for the quantum disguise 
trapping potential, the optimization procedure has a ten- 
dency to converge to a trivial solution, with the optimized 
parameters Xi diverging from the molecule position xq. 
In this case, if there is a negligible probability of electrons 
occupying Vmoh the presence or absence of the molecule 
does not change the response of the nanostructure. To 
avoid such solutions, we introduce a renormalized con- 
trast ratio Icontr/5p% where 5pg = Jy \pi{r)- p2{r)\'^d^r, 
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where pi and p2 are the electron ground state density 
with and without a molecule, respectively. 

In Figs. [2] (a)-(c) we present the calculations in the 
same order, as in Figs. [T] (a)-(c), but for the "quantum 
disguise" potential. For the "quantum disguise" problem 
we set the parameters of the molecular potential Vmoi (r) 
io A ^ 0.06e/i and B — 0.14i~^, so the perturbation 
due to the presence of the molecule becomes comparable 
to the original trapping potential Vtrap- Note, while the 
trapping potentials with and without the molecule, and 
the corresponding electron densities shown in Fig[21[c) 
are rather different, the optical responses of both systems 
at frequency lu = O.lE'o is very close to each other (see 
FigHd)). 

The effect of quantum disguise can be explained in 
terms of quantum interference between transitions in the 
combined system + Vmoi, which cancel each other 
to give approximately the response of 7J„. We checked 
this assumption by plotting the induced charge densities 
Pij (x) , corresponding to the transitions between states i 
and j. We found that all these densities are quite dif- 
ferent for the potentials shown in Figl^^a) and Figl^jjc). 
However, they all sum up to almost identical total in- 
duced charge densities Pind{x) = X^ij Piji^)- 

We also studied the dependence of the nanostructure 
response on the frequency of the external field. Fig. [3l^a) 
shows the contrast ratio as a function of the frequency of 
the incoming field lo. The dashed line shows Icontri^) for 
the initial, non-optimized trapping potential in Fig. [ija), 
and the solid line corresponds to the optimized trapping 
potential in Fig. \V[c). Note, that the solid curve has 
a maximum at the frequency hio — O.lii^O; a-t which the 
optimization was performed. In Fig. [3fb) we plot the 
contrast ratio Icontri'jj) for the potentials shown in Figs. 
[2ja),(c). Again, the dashed hue corresponds to the ini- 
tial, non-optimized trapping potential, and the solid line 
corresponds to the optimized one. Note the dip near 
the target frequency hco = O.IEq^ at which the optimiza- 
tion was performed. In both figures the improvement of 
the target functionality Icontri^) is of two-three orders 
of magnitude. 

We now consider a molecule position sensitivity analy- 
sis. In Figs|4l^a,b) we show the contrast ratio Icontr as a 
function of the molecule position for the non-optimized 
trapping potential and the optimized one. In FiglDJa) we 
show the results for the trapping potential parametrized 
using Np = 3, and in FiglD^b) Np = 5. It is clear, that 
even without any localization mechanism an array con- 
sisting of the optimized nanostructures will outperform 
a similar array of non-optimized ones. Indeed, assuming 
that the molecule's position xq is uniformly distributed 
on [0, L], one can estimate the avarage enhancement by 
integrating the area under the curve. The integration 
over the curves shown in FigHJa) gives the average value 
of the contrast ratio oi S = 7.01 with optimizion, and 
5" = 0.46 for a nonoptimized nanostructure. Integra- 
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FIG. 1: (a)-(b): Initial non-optimized trapping potential, (a) 
Trapping potential, no molecule (thin solid line) and with 
molecule (red dotted line). The corresponding electron den- 
sity, no molecule (thick solid line) and with a molecule (red 
dashed line), (b) Induced field intensity, no molecule, (solid 
line), and with a molecule (dotted line), (c)-(d) Same, as 
(a)-(b), but for an optimized potential that maximizes the 
contrast ratio Icontr at the frequency of the external field 
huj = O.l-Eo. Note for (d) we used a logarithmic scale for 
the j/-axis. 

tion of the curves shown in Figl4l(b), gives S = 8.43 and 
S' = 0.053 correspondingly. In the case of molecules 
localized near xq — L/2 the enhancement factor could 
reach up to three orders of magnitude, compared to the 
non-optimized nanostructures. 

In conclusion, we demonstrated the feasibility of op- 
timizing of the trapping potential parameterized by just 
a few control parameters. Even in this limited case the 
target functionalities (such as detector sensitivity) can 
be enhanced by orders of magnitude. 

Since the effects we consider are based on the resonance 
properties of the nanostructure, larger values of the level 
broadening 7 will result in reduced improvement of the 
design. From a simple analysis one can estimate that the 
intensity enhancement for the optimized design scales as 
oc (7/i?o)~^. Thus, phonon, electron-electron scattering 
and other processes that lead to level broadening will 
reduce the optimization effect. 

The results of optimal design for the quantum disguise 
problem show that although the real and imaginary parts 
of the nanostructure response are coupled through the 
Kramers-Kronig relations and are not completely inde- 
pendent, we succeeded to make the scattering and ab- 
sorbtion properties of two significantly different quantum 
systems almost identical at a given frequency . 

It is known that for the electron Hamiltonians there 
is a one-to-one correspondence (within an additive con- 
stant) between the non-degenerate ground state electron 
density p(r) and the trapping potential Vt,ap(r). This 
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dependence is established through the Hohenberg-Kohn 
mapping 1^ . Although some physical properties of the 
system can be deduced from its ground state, many char- 
acteristics, such as transport properties and optical re- 
sponse, strongly depend on the excitation spectrum of 
the system. The results suggest that there are situ- 
ations, when the precision of the calculated excitation 
spectrum is crucial. In the case of the single molecule 
detector shown in FiglU two quantum systems with very 
close ground state densities could have very different re- 
sponses. And the results for the quantum disguise prob- 
lem show that two quantum systems with almost iden- 
tical optical response at a given frequency could have 
rather different electron ground states. 



FIG. 2: Initial non-optimized trapping potential. (a) 
Trapping potential, no molecule (thin solid line) and with 
molecule (dotted line). The corresponding electron density, 
no molecule (thick solid line) and with a molecule (dotted 
line), (b) Induced field intensity, no molecule, (solid line), 
and with a molecule (dotted line), (c)-(d) Same, as (a)-(b), 
but for an optimized potential that minimizes the contrast 
ratio Icontr at the frequency of the external field huj = O.IEq. 



The relatively high sensitivity of the optical response 
of the optimized nanostructures to fabrication errors can 
potentially pose a serious problem. To address this prob- 
lem we suggest that the optimal properties of the nanos- 
tructures could be controlled during the assembling phase 
using closed loop adaptive techniques similar to those 



used in optical field control [13 1. 
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FIG. 3: The contrast ratio Icontr between the induced field 
intensity, with and without the molecule, as a function of the 
frequency of the external field ui. (a) The dashed line refers to 
the non-optimized trapping potential shown in Fig. [2a), and 
the solid line is for the optimized potential in Fig. [He). Note 
the local minimum near the target frequency hu — 0.1i?o- 
(b) Same as (a), but the contrast ratio is calculated for the 
potentials shown in Fig. [U^a) and (c), correspondingly. 




FIG. 4: The contrast ratio Icontr between the induced field 
intensity, with and without the molecule, as a function of 
position xo of the molecule, (a) The dashed line refers to the 
non-optimized trapping potential shown in Fig. [ija) , and the 
solid line is for the optimized potential in Fig. [ijc) . (b) Same 
as (a), but the contrast ratio is calculated for an optimized 
potential using {xi}, i — 1,..,5 parameters. 
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